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Abstract 

Complete discrimination system for polynomial and direct integral 
method were discussed systematically. In particularly, we pointed out 
some mistaken viewpoints. Combining with complete discrimination sys- 
tem for polynomial, direct integral method was developed to become a 
powerful method and was applied to a lot of nonlinear mathematical 
physics equations. All single traveling wave solutions to theses equations 
can be obtained. As examples, we gave all traveling wave solutions to 
some equations such as mKdV equation, Sine-Gordon equation , Double 
Sine-Gordon equation, triple Sine-Gordon equation, Fujimoto-Watanabe 
equation, coupled Harry-Dym equation and coupled KdV equation and 
so on. At the end, we give the partial answers to Fan's a problem. 
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1 Introduction 

In order to solve the traveling wave solutions to some nonlinear evolution equa- 
tions, sometimes these equations are reduced to an ordinary differential equa- 
tions such as 

u'(t) = G(u,0 1 ,--> ,0 m ), (1) 
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where 0\, ■ ■ ■ ,8 m are parameters, its solutions can be wrote by integral form 



Thus what we need is only to solve the integral (2). According to different pa- 
rameters, we will give different solution to the integral. This is so called direct 
integral method. Although direct integral method is a routine method, but to 
decide the parameter's scope is still rather difficult. Thus the most important 
steps are to decide the parameter's scopes and the corresponding integral so- 
lutions. A new mathematical tool named complete discrimination system for 
polynomial [1-5]) is applied to this problem so that the parameter's scopes is 
solved. 

It is just because of combining with complete discrimination system that 
direct integral method becomes a powerful and efficient method. There are 
a lot of nonlinear mathematical physics equations such as KdV equation, Hi- 
rota equation, breaking soliton equation, Coupled mKdV equation, Chen-lee-Liu 
equation, Kundu equation equation, cubic Schrodinger equation, Sine-Gordon 
equation, double Sine-Gordon equation, Sinh-Gordon equation, double Sinh- 
Gordon equation, triple Sinh-Gordon equation, Combine KdV-mKdV equation 
and nonlinear Klein-Gordon equation and so on, their all traveling wave solu- 
tions can be obtained by using direct integral method. These complete results 
can't be obtained by other any indirect methods ( see [6-19]) such as trial equa- 
tion method([19]), Jacobi elliptic function expanse method ([12]), sub-equation 
method([15,16] and homogenous balance method and so on. 

Although direct integral method is the most basic method, but there are 
two obvious mathematical facts about the integral(2) to be ignored. One fact is 
that different variable transformation don't offer to a different solution, different 
variable transformation only give different representation of solution. Another 
fact is that it is just different parameter's scope to offer to different solution. 
Some authors didn't notice the fact onc( [20-23]), they still think that a new 
variable transformation can bring new solutions. We must point out that this 
is a mistaken viewpoint. In fact, the elementary calculus tell us that a defined 
function has infinite primitive functions which are different each other by a con- 
stant, and that a new variable transformation only offer to a new representation 
of the primitive function. 

In the present paper, using complete discrimination system for polynomial 
and direct integral method ,we give the classifications of all single traveling wave 
solutions to a lot of nonlinear differential equations such as mKdV equation 
and Sine-Gordon equation, Fujimoto-Watanabe equation, coupled Harry-Dym 
equation and coupled KdV equation and so on. An the end, we give the partial 
answer to Fan's a problem. 

This paper is organized as follows. In section 2, we deal with some equations 
whose reduced ODE's take the form of (u 1 ) 2 = a 3 u 3 + a 2 u 2 + a\u + a , and ob- 
tain the classifications of all single traveling wave solutions to those equations. 
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In section 3, we deal with some equations whose reduced ODE's take the form 
of (u 1 ) 2 = a 4 u 4 + a^u 3 + a 2 u 2 + a\u + ao, and obtain the classifications of all 
single traveling wave solutions to those equations. In section 4, we give some 
applications of the complete discrimination system for the fifth order polyno- 
mial. In section 5, we deal with two equations with rational form of reduced 
equations and obtain classifications of their all single traveling wave solutions. 
In section 6, we use a one-parameter's thick obtain the classifications of all single 
traveling wave solutions to some equations such as coupled KdV equation and 
so on. In section 7, we discuss Fan's two problems, and give an answer to the 
first problem. 

2 Applications of the complete discrimination 
system for the third order polynomial 

There are a lot of nonlinear evolution equations such as KdV equation 

u t + 6uu x + u xxx = 0; (3) 

BBM equation 

u t + u x + uu x + au xxt = 0; (4) 

Klein- Gordon equation 

u u ~ u xx + au + f3u 3 = 0; (5) 
Landau-Ginburg-Higgs eqution 

u t t - u xx - m 2 u + g 2 u 3 = 0; (6) 

Sine-Gordon eqution 

u xt = asinu; (7) 
Dodd-Bullough-Mikhailov equation 

u xt + v exp(w) + q exp(-2w) = 0; (8) 

Zhiber-Shabaut equation 

Uxt + P cxp(w) + q exp(-u) + r exp(-2u) = 0; (9) 

Sinh-Gordon eqution 

u xt = asinhu; (10) 
Joseph-Egri equation Sine-Gordon eqution 

u t + u x + uu x + 0u x tt = 0; (11) 

Kundu equation 

m - u xx - 2i(2a - l)\u\ 2 u x * +a(Aa - l)|u| 4 w = 0; (12) 
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derivation Shrodinger cqution 

u t = m xx + i(c 3 |u| 2 + c 5 |u| 4 )w + s 2 {\u\ 2 u) x = 0; (13) 
Ablowitz equation 

m = Uxx - ii{\u\ 2 u x * +8|u| 4 u; (14) 
and so on. They all can be reduced to the following third order integrable ODE 

(u 1 ) 2 — a 3 u 3 + a 2 u 2 + a,\u + ao, (15) 

So if we give the classification of all solutions of the Eq.(15), then all single 
traveling wave solutions to those equations (3)-(14) can be classified. According 
to the Ref.^1) we list those results for convenience in the following, let 

w=(a 3 )^u, d 2 = a 2 (a 3 )~i , di = ai(a 3 ) - s, d = a . (16) 

Then Ea. (|2T3|) becomes 

±(a 3 )%-£o)= f . 3 I dw. (17) 

J Vw J + d 2 w z + d\w + do 

Denote 

F(w) = w 3 + d 2 w 2 + d x w + d - (18) 
A = -27(ff + do-^) 2 -A(d 1 -§) 3 , (19) 

Di = di - |, (20) 

where A and Di make up a complete discrimination system for F(w). There 
are the following four cases to be discussed(jTHj): 

Case 1: A = 0, D 1 < 0, then we have F(w) = (to - a) 2 (uj - /?),a ^ (3. If 
w > f3, the solutions are as follows: 



U =(a3)-^(a-/3)tanh 2 (^^( a3 )^(e-eo))+/3], a > (3; (21) 



U =( a 3)-^(«-/3)coth 2 (^^(a 3 )^(C-eo)) + /3], «>/3; (22) 
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« = (as)"* [(-« + 0) sec 2 (^4±^(a 3 )^ (£ - &)) + /3], a < (3. (23) 
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Case 2: A = 0, D\ = 0, then we have F(w) = (w — a) 3 , the solution is as 
follows: 

U = 4(a 3 )-§(£-£o)- 2 + a. (24) 



Case 3: A > 0, Di < 0, then F(w) = (w — a)(w — fi){w — 7), we suppose 
that a < (3 < 7. When a < w < (3, we have 



u = (a 3 )-*[a + 09 - a) sn 2 (^-^(a 3 )^ - £o)),m)]. (25) 



When w > 7, we have 

u = (a 3 )~3 

where m 2 = ^1^. 

7 — a 



, 7 -/?sn 2 (^£(a 3 )^-eo)),m) 
cn 2 (^(a 3 )^(e-eo)),m) 



(26) 



Case 4: A < 0, then we have F(w) = (w — a)(w 2 + pw + q), p 2 — 4q < 0. 
We have 



/ r^> 2W~a*+pa + q 

= (a 3 ) 3[a-v a +P a + ^ /^ 77— r 



2 1 r 1 1 ], (27) 

1 + cn((a 2 +pa + ?) 4 (a 3 )3 (^ - ^o),m) 

where m 2 = i(l , " +2 ). 

We take only three equations as examples to illustrate this method, other 
equations can be dealt with similarly. 

Example 1. Classification of all single traveling wave solutions to KdV equa- 
tion 



Although KdV equation is a classical equation, its single traveling wave 
solutions had been discussed many years ago. But its classification hasn't been 
reported, some author also try to give its new single traveling wave solutions, 
so we take it as an example and give its all single traveling wave solutions. In 
fact under the traveling wave transformation, KdV equation can be reduced to 
the following ODE: 

( W ') 2 = -^ u3 -^3" u2 + ClM + C "- ( 28 ) 

Thus it is easy to give the classification of all solutions to ODE 1281 according to 
the former results. For briefly we omit them. 
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Example 2. Classification of all single traveling wave solutions to Sine- 
Gordon equation. 



Sine-Gordon equation reads 

u xt = asinu, (29) 

which widely applied in physics and engineering. Its so-called new traveling 
wave solutions have been obtained in different functional forms by different 
methods([25-30]). Among these some method([30]) is very complex but not 
essence and natural. But we must point out that it is very easy to give all 
traveling wave solutions to Sine-Gordon equation. In the following, we only 
need elementary integral method to do this thing. Instituting traveling wave 
transformation u = w(£), £ = kx + cot into Ea. (|29|l and integrating once yield 

^ / V—tiwau-c)' (30) 
where c is an arbitrary constant. Take variable transformation 

u = arccosw, (31) 

then Ea. (|30|l becomes 

± ^ " = / fa, w r ( 32 ) 

J ^-^(w-l){w + l){w-c) 

Firstly we assume -j^ > 0. There are the following three cases to be discussed: 
Case 1: c = 1, then we have 

u = arcsm{2tanh 2 (iw-^(£- £ )) - 1}; (33) 

,1 [2a 
3th' ( 

Case 2: c = — 1, then we have 



u = arcsin{2cotli (-W — (£ - &)) - 1}- (34) 



u = arcsin{-2sec 2 (-!-W-^(£ - f )) + !}• (35) 
2 V few 

Case 3: c ±1, we reorder 1, —1, c and denote them by a < (3 < 7. In fact 
there are three case such as 1 > c > — 1, c > 1, and — 1 > c. When a < w < (3, 
we have 



u = arcsin{a + (J3 - a) sn 2 (^A/ ^(£ - £„)),m)}. (36) 



2 V few 
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When w > 7, we have 



.7-/3sn 2 (^Jf^-£ )),m). 



u = arcsin{ yku, ^ 

™ 2 (^\/iH£-£o)),m) 



where m 2 = — — . 

*y — ct 



If r 5 - < 0, we take w — —v, then Ea.i|52)l becomes 

V a/-£t (v- 1 w + 1 w + c 



We have similar three cases: 
Case 1: c = 1, then we have 



u = - arcsin{2tanh 2 (i^--^(£ - £„)) - 1}; (39) 



= - arcsin{2coth 2 (i^-^(£ - Co)) - !}• (40) 



Case 2: c = — 1, then we have 



u= -arcsin{-2sec 2 (iy--^(e-eo)) + l}. (41) 

Case 3: c 7^ ±1, we reorder 1, —1, c and denote them by a < (i < 7. When 
a < w < (3, we have 



2/V7- " / 2a , 



arcsin{a + (/3 - a) sn 2 (^ - — (£ - &)),"*)}. (42) 



When w > 7, we have 



7 _^ sn2( ^^: K _ ?o)) , m) 

it = — arcsmj }, (43) 

where m 2 = 

7 — a 

Example 3. Classification of all single envelope traveling wave solutions to 
Hirota equation 

Hirota equation reads 

iu t + u xx + 2\u\ 2 u + iauxxx + 6ia\u\ 2 u x = 0, (44) 
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when a = 0, it becomes Schrodinger equation. Take envelope travelling wave 
transformation u = u(£) exp(^), £ = x + cot, rj = px + qt, then separating the real 
part and the imagine part of Eq. I|44|> yields 

{uj + 2p - 3ap 2 )u' + 6au 2 u + au'" = 0, (45) 

(ap 3 -q- p 2 )u + 2(1 - 3ap)u 2 + (1 - 3ap)u" = 0. (46) 

In order to make the Eq. l|45[) and the Eg. (|46fl consistence, the parameters must 
satisfy the following condition: 



uj + 2p - Zap 2 



ap 3 — q — p 2 1 — Sup 
Under the condition (|47H . the above two equations become one equation 

/ l\2 4 . OLp 3 q p 2 2 

(u ) = -u H u + c 

3ap — 1 

where c is an arbitrary constant. The corresponding integral form is 



(47) 



±K - fo) = / ! f I ( 49 ) 

J . /_„4 I «P 3 -g-p-„,2 I r . 



-« 4 + ^k^ u/ + c 

in order to solve the above integral, we take the change of variable w = u 2 , then 
the integral becomes 

±2(C - Co) = / , ^ a (50) 

Therefore according to the classification of all solutions of the Eq.(15), it is easy 
to give classification of all single envelope traveling wave solutions of the Hirota 
equation. We omit concrete expressions for briefly. 

3 Applications of the complete discrimination 
system for the fourth order polynomial 

We first write the complete discrimination system for the fourth order polyno- 
mial f(w) = w A + pw 2 + qw + r as follows(0]): 

Di = 4, D 2 = -p, D 3 = 8rp - 2p 3 - 9q 2 , 

27 

Di = 4/r - p 3 q 2 + 36prq 2 - 32r 2 p 2 - —q 4 + 64r 3 , 

F 2 = 9p 2 - 32pr. (51) 
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Then we consider the following ODE 

(w'(0) 2 ^ e ( w4 +PW 2 + qw + r), (52) 
where e = ±1. Rewrite Ea. (|52[) by integral form as follows: 

±(£ - &) = / dw (53) 

J e(w 4 + pw 2 + qw + r) 

Although Eq.(52) has been studied extensively, some solutions were obtained 
under some special parameter's conditions. But we will give in the following 
the classification of its all solutions. According to the complete discrimination 
system for polynomial F(w) = w A + pw 2 + qw + r, there are nine cases to be 
discussed ( "4), we list those results for convenience in the following. 

Case 1: L> 4 = 0, D 3 = 0, D 2 < 0. Then we have 

F(w) = ((w-h) 2 + s 2 ) 2 , (54) 
where h, s± are real numbers, si > 0. When e = 1, we have 

to = s 1 tan[si(£-£ )]+Zx- (55) 



Case 2: D 4 = 0, D 3 = 0, D 2 = 0. Then we have 

F(w) = w 4 . (56) 

When e = 1, we have 

1 



(57) 



Case 3: D 4 = 0, D 3 = 0, D 2 >0,E 2 = 0. Then we have 

F(w) = {w-a) 2 (w- P) 2 , (58) 
where a, (3 are real numbers, a > (3. If e = 1, when w > a or w < [3, we have 

^ = ^[coth^(£-£o) -l]+/3. (59) 
When a > w > (3, we have 

™=^[tanh^(£ -&)-!]+/?. (60) 
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Case 4: D A = 0, D 3 > 0, D 2 > 0. Then we have 

F(w) = (w-a) 2 (w-(3){w-- f ), (61) 

where a, (3, 7 are real numbers, and (3 > 7. If e = 1, when a > /? and w > 8, or 
when a < 7 and w < 7, we have 

± ( _ } = 1 ^EEEWE^E^EWEEL. (62) 

y/(a-0)(a-j) \w-a\ 
When a > 8 and w < 7, or when a < 7 and w < 8, we have 

± ( _ } = 1 ^EEEME^EEEEEEEl. (63) 

V(« -/?)(<* -7) \w-<*\ 
When (3 > a > 7, we have 

, /, , x 1 • (w - (3)(a - 7) + (a - p)(w - 7) ,„„, 

± (g - go) = ; arcsin ^ ^-7 LZ. (64) 

V(/3-a)(a- 7 ) -<*)(/?- 7 )| 



If e = —1, when a > /? and tu > /3, or when a < 7 and to < 7, we have 

± ( _ } = 1 ^VEEBEEE^EWEEl. (65) 

^(a -/?)(<* -7) 
When a > 8 and w < 7, or when a < 7 and w < /3, we have 

± ( _ } = i ^EEEMEE^EWEEEEEl. (66) 

^(a -/?)(<* -7) 
When /? > a > 7, we have 

, y s 1 • (w-/3)(a-7) + (Q!-/3)(w-7) 

± (g - Co) = = arcsin ^ ^-7 ^ ^. (67) 

V(/3-a)(a- 7 ) | -«)(/?- T ) | 



Case 5: D 4 = 0, D 3 = 0, D 2 > 0, £ 2 = 0. Then we have 

F{w) = (w-a) 3 (w- 0), (68) 

where a, (3 are real numbers. If e = 1, when w > a, w > (3, or w < a, w < (3, we 
have 

W= (a-mg-6.) a -4 - (69) 
If e = — 1, when w > a, w < (3; or w < a, w > (3, we have 

4(8 -a) , . 

W= 4 +(Q - m C-^ - (?0) 
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Case 6: D 4 = 0, D 2 D 3 < 0. Then we have 

F(w) = {w-a) 2 {{w-h) 2 + s\), 
where a,l\ and Si are real numbers. If e = 1, we have 



w = 



where 



exp [±y/(^WT4(t - &)] -j+y/(a- h) 2 + s( 

a - 2Zi 



7 : 



^( a _; l)2 + s 2- 



Case 7: D 4 > 0, D 3 > 0, L>i > 0. Then we have 

F(w) = (w — ai)(w — a 2 )(w — a 3 ){w — a 4 ), 

where a\, a 2 , a 3 , a 4 are real numbers, and ot\ > a 2 > a 3 > a 4 . If e = 
w > a.\ or w < a 4 , we have 



w 



a 2 (a 1 -a 4 )sn^(-^ ^ -{€-£o),m 



(ai - a 4 )sn 2 3 (£ _ £o),m 



When a2 > w > «3, we have 



a 4 (a 2 - a 3 ) sn^ (-^ ^ -(£ - £o),m 



(a2 - 03) sn^ ' 



2 "(£-&>), m 

If e = — 1, when ai > w > a2, we have 



w = 



a 3 («i ~ « 2 ) sn 2 ( V^E^E^ _ £ o) , 



/ \ 9 / \l (ai-oa)(a2-a4) />. >. n 

(ai - a 2 )sn 2 ^ i(£ - £ ), 

When a 3 > w > a 4 , we have 

/ x 2 / V (ai-a 3 )(a 2 -a 4 ) , f £ , 

t(a 3 - a 4 )sn z j (s - So), 



w = 



rn 



I \ 9 / v/ («l-«3)(a2-« 4 ) //■ f \ 

(a 3 - a 4 ) sn 2 (-^ ^ H£ - £ ), m 



— a 4 (a2 — a 4 ) 



- (a 2 - a 4 ) 



a 3 (a:2 — a 4 ) 



(a 2 - ck4) 



- a 2 (a x - a 3 ) 



(ai - a 3 ) 



a 4 (a 3 - ai) 



(a 3 - ai) 



where m 2 



(ai-tta)(tt3-tt4) 
(«l- a 3 )(a 2 — out) ' 



Case 8: L> 4 < 0, D 2 D 3 > 0. Then we have 

F(w) = (w- a){w - p)((w - h) 2 + s\), 
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where a,(3,h and si are real numbers, and a > [3 >,s± > 0, we have (here 
positive sign for e = — 1, negative for e = 1.) 



ccn ( 2mmi (f-£o),m)+d 



(80) 



2mmi 

where 

c = a — l\ — , d = a — l\ — s\mi, 

mi 

a=±[(a + p> -(a- (3)d],b = i[(a + 0)d -(a- p>], 

£? = ^ + (a-w-/i) mi= £ ±v /^n v = T J-3, (81) 

si(a — p) 1 + mf 

we choose mi such that em! < 0. 

Case 9: D4 > 0, D2D3 < 0, then we have 

F(w) = ((w - hf + 4)(w - l 2 ) 2 + si), (82) 
where Zi, Z2, s\ and s 2 are real numbers, and s x > s 2 > 0. If e = 1, we have 

^ = asn(?7(g~g ),m) + &cn (?7(g-go),m) 
csn ~ Co), m) + den (??(£ - Co), m) ' 



where 



s 2 

c = — si , d — h—l2, a = l\c + s\d, b = hd—s\c, 

mi 

E = Z + Z + ^-V\ mi =E + y/&=l, 



1 m 2 c 2 + d 2 



Remark 1: For case 7 and case 8, we can use another transformation to 
obtain more simple representations. For example, for case 7, we consider only 
e = 1, as an example, we take a transformation as follows: 

wi = , (85) 

w — Ct\ 

where A — (ai — Q2)(o!i — a.z){a.\ — 0:4), then corresponding integral becomes 

1 , n f dwi 
±a*(Z-io) = — 86 
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where 0i = A? , 3 2 = A ^ , 3 3 = A? . when A > Wi > 8 2 , we have 



= /9a + (/3 2 - /3 3 )sn 2 ( v ^ " J (£ - &),m). (87) 



When toi > /?i , we have 



/3i-/3 2 sn(^3I(£-£o),m) 



= - ^ 2 " — , (88) 
where m = o 2 ~f! 3 . 

Px—ps 

For the case of e = — 1 and the case 8, we can deal with them similarly. 

Remark 2: All above integrals can be fond in Ref. 24 . 

Example 1. Classification of all single travelling wave solutions to mKdV 
equation. 

mKdV equation reads 

U t + au 2 u x + bu xxx = 0. (89) 

Its exact traveling wave solutions have been studied extensively ([20-23]). But its 
all traveling wave solutions haven't been reported. Using direct integral method 
and complete discrimination system for the fourth order polynomial, we give in 
the following all traveling wave solutions to mKdV equation. 
Take traveling wave transformation as follows: 

u = u(£),£ = kx + wt. (90) 

Instituting Eq. (|90|l into Ea. H89l) and integrating twice yield 

{u'f = A(u 4 +pu 2 +qu + r), (91) 

where 

A a a ^ (92) 

and q, r are two arbitrary constants. Write Ea. (|91|l by integral form 

±V^i(£ - &) = / dU (93) 

where e = ±1. If A > 0, we take e = 1. If A < 0, we take e = — 1. According 
to the former results, we can give the classification of all single traveling wave 
solutions to MKdV equation, we omit them for briefly. 
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Example 2. Classification of all single traveling solutions to combined 
KdV-MKdV equation([m,[H ]) which reads 

u t + auu x + bu 2 u x + u xxx = 0, (94) 

where a and b are physical parameters. Taking traveling wave transformation 
u = = kx — tot, and instituting it into Eq.JMJ) and integrating once 

yields: 

(u 1 ) 2 = a^u + a 3 u 3 + ci2U 2 + a±u + ao, (95) 

where 

°* = -fl^. a ^'W a2 = ¥> (96) 
and Oi, ag are arbitrary constants. If a < 0, we take transformations as follows: 

w = al(u+^-), £ = aiti. (97) 

Then Ea. (|95(l becomes 

{w') 2 =w A + pw 2 + qu + r, (98) 

where 



a 2 -\. a\ a 2 a 3 aia 3 a 2 a% 3af 

y/ol oa^ 2a4 4a4 loa^ 25604 



If a > 0, we take transformations as follows: 



w= {-a A )*(u+ -?-), £=(-a 4 )^i 
4a4 



(100) 



Then Ea. (|93|l becomes 



(w') 2 — ~(w 4 + pw 2 + qu + r), 



(101) 



where 



P= q=(-a4) 4 (-^r + ^ ai). r = -a - 

8a| 2a 4 



-a4 



0103 a2 a i 3a| 



4a 4 16a| 256a 4 i ' 
(102) 



According to the classification of solutions to ODE (w 1 ) 2 = ±(w 4 + pw 2 + 
qw + r) , we can give all traveling wave solutions to combined KdV-MKdV equa- 
tion. We omit them for simplicity 
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Example 3. Classification of all single traveling wave solutions to Konopelchcnko- 
Dubrovsky equation which reads 

3 

Ut = u xxx - 6f3uu x + -a 2 u 2 u x - 3w y + 3auu x = 0, (103) 

w x = u y . (104) 

Takeing travelling wave transformation u = u(£),£ = k\x + k^y + u)t, institut- 
ing it into Ea. (|103[) and integrating twice yield the same form equation with 
Eq.JHUl: 

(u 1 ) 2 — a^u A + a^u 3 + ci2u 2 + a\u + ao, (105) 

where 

a 2 a -2/3 2(k 1 u; - 3k 2 ) 
fl4 "4fc2' a3 = ^fcf~' a2 " kf ' ( ' 

and ai,ao are arbitrary constants. Notice > 0. According to the example 1, 
it is easy to write all traveling wave solutions to KD equation, we omit them 
for briefly. All nine cases under the condition 04 > in section 4 are possible. 

Example 4: Classification of all single traveling solutions to double Sine- 
Gordon equation f |42U43| ) which reads 

u xx — u t t — sin au + A sin — u, (107) 

where a and A are nonzero constants. Its some exact traveling wave solutions 
have been obtained. We give in the following its all traveling wave solutions. 

Take traveling wave transformation 

u = u(£ 1 ), £ 1 = kx-0Jt, k^±u. (108) 



Instituting Eq. (|108|) into Eq. l|107fl and integrating once yields 

±(fi " £10) = / , ^ (109) 

J ^Jc2 — c\ cos au — Co cos -^u 

where 

2 4A 
Cl = ~Ti 2\ ' c ° = ~Ti 2V ( U °) 

In order to solve above integral (|109fl . we take transformation as follows: 

2 , . 

u = — arccosu. (HI) 
a 
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Instituting Eq. (|lll|) into integral (|109|> . we have 

f Idv 



If ci > 0, let 



y/{v -1)(V+ l)(2 Cl V 2 + C Q V + Cl - C 2 ) 



I— i c ° c a c 
4ci V 2c i 



Thus Ea. l|112|l becomes respectively 

±(e-&) = / 

where 



^/(w; — a)(w — (3)(w 2 + d ) 



c c cx-c 2 c 



If a < 0, let 



i c ° t a c 

w = y—civ — - — , £ = — =£1. 

4ci a/2ci 



Thus Ea. (jll2f) becomes respectively 

±(£-G>) = / 

where 



cIj; 



Co . > „ c , ci - c 2 c 



Denote 



f(to) = ( w - a )(u> - /3)(w 2 + d ), 



and rewrite it as follows: 



where 



F(w) — w A + pu> 2 + gu> + r, 



p = do — (a+fl), r = —do(a + /3), s — d ap. 



But because of special form of F(w) (see Ea. ()119(l ^. we know that the 
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case 2, case 6 and case 9 in section 3 are impossible to occur in the case of 
double Sine-Gordon equation. Other cases are possible. Thus we can give all 
traveling solutions to double Sine-Gordon equation according to the cases 3-5 
and cases 7-8 in section 3. We omit them for briefly. 

Example 5. Classification of all single traveling wave solutions to Pochharmer- 
Chree Equation which reads 

u tt - au xx + f3u + 7U 3 + Su 5 . (122) 

Takeing traveling wave transformation u — u(£), £ = kx + u>t, instituting it into 
Ea. l|122fl and integrating once yields 

(u r ) 2 = a 3 u 6 + a 2 u 4 + aiu 2 + a , (123) 

where 

a 3 — ~TT~5 TIT' ° 2 ~ _ o7^ TIT' ai ~ 2 uz> \ Ll V 

and ao is an arbitrary constant. Let v — u 2 , then Ea. (|123| l becomes 

{v'f = 2v(a 3 v 3 + a 2 v 2 + a x v + a ). (125) 

Similar to the example 3, we can give all traveling wave solutions to Pochhamer- 
Chree equation (|122fl . Because of special form of the right polynomial of 
Ea. l|125|l . the case 1 and case 9 in section 3 are impossible to occur, all other 
cases 2-8 are possible. We omit them for simplicity. 

Example 6: Classification of all single traveling wave solutions to classical 
Bousinesq system 

Vt + [(1 + V)u]x = -^u xxx , (126) 



u t + uu x r) x = 0. (127) 

Take traveling wave transformation 77 = ??(£), u = u(£), ^ = kx + tot, instituting 
it into above system and integrating once yield 

fc 3 

c +w77 + fc(l + ?7)M = -— u", (128) 

V = Cl ~yu-lu 2 . (129) 
fc 2 

Instituting Eq. (|129|) into Eq. l|128|l and integrating once,we have 

(u') 2 = a4« 4 + a 3 u 3 + (I2U 2 + a\u + ao, (130) 
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where 

1 8w 4^ 2 -4fc 2 (l + Cl ) 
fl4= p' a3 = 5fc5' ° 2 = £3 ' ai=c +ci^, (131) 

and ao is an arbitrary constant. Notice en > and according to the example 
1, we can easily to obtain all traveling wave solutions to classical Boussinesq 
system. For briefly we omit them. 



4 Applications of the complete discrimination 
system for the fifth order polynomial 

We consider the following ODE 

(w'(0) 2 =F(w) = w 5 + pw 3 + qw 2 + rw + s. (132) 

We write its complete discrimination system as follows (see PQE]): 

D 2 = -p, D 3 = 40rp - 12p 3 - 45q 2 , 
D 4 = 12p 4 r - 4p 3 q 2 + \Ylprq 2 - 88r 2 p 2 
-AOqsp 2 - 27q* - 300qrs + 160r 3 , 
D 5 = -WOOqsr 3 - 3750pqs 3 + 2000ps 2 r 2 - Ap 3 q 2 r 2 + 16p 3 q 3 s 
-900rsV + 825pVs 2 + U4pq 2 r 3 + 2250rq 2 s 2 + 16p 4 r 3 
+108p 5 s 2 - 128rV - 21r 2 q 4 + 108sg 5 + 256r 5 + 3125s 4 
—72rsqp i + 560sqr 2 p 2 — G30prsq , 

E 2 = 160rV + 900gV - 48rp 5 + 60rP 2 q 2 + 1500pqrs + 16<?V 

-llOOqsp 3 + 625s V - 3375sg 3 , 

F 2 = 3<? 2 - 8rp. (133) 

According to the complete discrimination system (|133|) for polynomial f(w), 
there are twelve cases to be discussed. 

Case 1: D 5 = 0, D 4 = 0, D 3 >0,E 2 ^ 0. Then we have 

F{w) = {w-a) 2 (w- 0) 2 {w--/), (134) 
where a, 8, 7 are real numbers, a 7^ (3 7^ 7. When w > 7, we have 



a — 3,. . , -v/w — 7 / -i/u; — 7 

± - (g - Co) = V7 " Q arctan - V7 - P arctan —===, 

2 V7 ~ " V 7 - P 

7>a, 7 >/3; (135) 



± — - — (4 — £ ) = v 7 — a arctan 
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xh| ^;^2 |, 7 >a, 7 <ft (136) 
V w 7 + V P - 7 



2 2V^7 V^ T+ y^ T 1 2v? _ 7 

xln| >fz2-^Z2 |, 7 <a, 7 </3. (138) 
^/to - 7 + VP - 7 



Case 2: L> 5 = 0, L> 4 = 0, £> 3 = 0, L> 2 7^ 0, F 2 7^ 0. Then we have 

F(w) = {w- af(w - /3) 2 , (139) 
where a, /3, are real numbers, a ^ [3. When it; > a, we have 



±— -K-&) = j— - V^-/3arctan^=,a> # (140) 

2 — a V a — P 



, a ~ ft a\ 1 1 . \[w~^a — \fft~~a. „ < _! . I \ 
± „ (£ - Co) = ; In - Vo. < /5. (141) 



Case 3: D 5 = 0, L> 4 = 0, £> 3 = 0, L> 2 ^ 0, F 2 = 0. Then we have 

F(w) = {w-a) 4 (w- p), (142) 
where a, p, are real numbers, a 7^ p\ When w >a, we have 

± —«"*>) = -2(^0) ^fa==p arCtan 77^' « < A ( 143 ) 



j_ a -ft(c c \ Vw-a 1 . /3 \/o; /? . . 

2 2(w — a) 4va — p vw — p + \Ja — p 



Case 4: L> 5 = 0, D 4 = 0, £> 3 = 0, D 2 = 0. Then we have 

F(w) = {w - a) 5 , (145) 
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where a is real number. When w > a, we have 



(146) 



Case 5: D 5 = 0, D 4 = 0, L> 3 < 0, £ 2 ^ 0. Then we have 

F(w) = (w-a)(w 2 +rw + s) 2 , (147) 
where a is real number, r 2 — 4s < 0. When w > a, we have 



± « " = " ^=7^ (cos * arctan w -.-p 2 + — 

W -«-p?-2pcoWji>]^ 

io — a — p 2 + 2pcos<pvw — a 



where 



1 v/4 



s — r 
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p = (a 2 +ra + s)*,(/) = - arctan ^ — . (149) 

2 —2a — r 



Case 6: D 5 = 0, D 4 > 0. Then we have 

F(w) = (to — a) 2 (w — ai)(to — 02) (to — 03), (150) 
where a, ai, 02 j «3 are real numbers, cti > a.2 > a 3 . We have 

±(£ - £ ) = -7 ^ {F(y, fc) - ^^1%, fc)}, (151) 

(a — a2jv a 2 — a 3 a\ — a a.\ — a 

where a ^ a\,a ^ a2,cc ^ Q-s and 

F(p,fc) = [ V ^ a , (152) 

1%, », fc) = f , 2 = , 2 ■ (153) 

Jo (1 + nsin ip) v 1 — k z sm <p 

Case 7: D 5 = 0, D 4 = 0, L> 3 < 0, E 2 = 0. Then we have 

F(w) = (w-afdw-hf + sj), (154) 
where a, l\ and s\ are real numbers. When w > a, if a ^ Zi + si, we have 
tan 6* + cot si tan# + si cot 
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" « fa ,o"^,'-a) 5 in ~ v V' - e * + F(„. t ) - E(„, t) ), (155) 



if a = Zi + si, we have 



/ sin 3 20 1 

±(£-£o) = W^^(^arcsin(/csin^) -F(^,fc)), (156) 



where 



tan2<9 = ^— , fc = sin0, < < ^, E(<p,k)= J 1 - fc 2 sin 2 ^d^. 
a — ti 2 Jo 

(157) 



Case 8: D 5 = 0, D 4 < 0. Then we have 

F(w) = (w- a) 2 (w - /?)((«, - Z x ) 2 + s 2 ), (158) 

where a,Zi and si are real numbers. When w > (3, if a ^ l\ — Sitan0, a ^ 
h + si cot 0, we have 

^_ tan + cot Si tan 6* + Si cot 

± ^ ~ ^ ~ 2(sit&n9-l-a)y/^= ^' ^ " si cot 6> + Z, +a 

X % Z sin , * ^ + n<P, k) - Efr, k)}; (159) 

if a = Zi — Si tan we have 



±(C-£o) = J^^(^arcsin(/csin^) -F(^,fc)); (160) 



if a = Zi + si cot 0, we have 



/sin 3 20 1 v 7 !-^ 2 sin 2 v5 + Vl - A: 2 sin^ 

- &) = V "^ (F fc) - 7t=P ln ^ ); 

(161) 

where 

tan26>= — ^— , k = sin0, < < (162) 
P — h 2 



Case 9: D 5 =0,D 4 = 0, L> 3 > 0, E 2 = 0. Then we have 

F(w) = (w-a) 3 (w- ^{w-j), (163) 
where a, j3 and 7 are real numbers. When w > a > (3 > we have 



±— — = , E(arcsinW , J 

2 ^/a — 7 p^T V tt_ 7 
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-J , y (164) 

In other cases such &sw>(3>a>-f and so on, we can give the corresponding 
solutions similarly. We omit them for simplicity. 

Case 10; D 5 > 0, D 4 > 0, D 3 >0,D 2 > 0, we have 

F(w) = (w — ai)(w — a 2 ){w — a 3 )(w — a±){w — a 5 ). (165) 

Then the corresponding solutions can be expressed by hyper-elliptic functions 
or hyper-elliptic integral such as 

±(£-£o)= / ^ (166) 

J y/(w — ai)(w — ot2)(w — as)(w — on){w — a^) 

Case 11: D 5 < 0, we have 

F(w) = {w- ai )(w - a 2 )(w - a 3 )((w - hf + s\). (167) 

Then the corresponding solutions can be expressed by hyper-elliptic functions 
or hyper-elliptic integral such as 

±(g-6))= / , ^ (168) 

J y/ (w - ai)(w - a 2 ){w - a 3 )((w - h) 2 + s\) 



Case 12: D 5 > 0A(£»4 < V D 3 < V D 2 < 0), where A means "and", V 
means "or". Then we have 

F(w) = {w- a){{w - hf + sl)){{w - hf + si). (169) 

Then the corresponding solutions can be expressed by hyper-elliptic functions 
or hyper-elliptic integral such as 



±(£ -&)=/- . ^= (170) 



dw 

y/(w - a)((w - hfTsl))((w - l 2 f + 4) ' 



Example 1: All traveling solutions to triple Sine-Gordon equation([3],[42-45]) 

112 2 
u xx -u tt = csinu+ -asin-u+ -6sin-w, (171) 

where a, b and c are nonzero constants. Its some exact traveling wave solutions 
have been obtained. In Ref.[3], we used the similar idea to deal with it, but we 
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didn't recognize that its all single traveling wave solutions can be obtained. We 
give in the following its all traveling wave solutions. 

Take traveling wave transformation 

fi = kx - ut, k ± ±cj. (172) 

Instituting Eq. (|172fl into Eq. l|171|l and integrating once yields 

±(fi " £10) = / 1 dU (173) 

J ■ I c 2 — 2ci cos u — 2ai cos — 2b x cos |u 



where 



01 = -: 2' 5l = l 2 ' Cl ~ I " 2 ' ( 174 ) 

K 2 — L0 Z K 2 — UJ fc 2 — U) 



In order to solve above integral (|173fl . we take transformation as follows: 

u = 3arccosv. (175) 

Instituting Eq. (|175|) into integral l|173|) . we have 



±(& " £10) = J 



y> - l)(v + l)(8ciw 3 + 4&iw 2 + (2ai - 6a)v - 2b x - c 2 ) 

(176) 



Let 

„ = ( 8ci )-*«,- A., C = |(8 Cl )*Ci- (177) 

Thus Ea. l|176|l becomes 

±(£-&)= / ; ^ . (178) 

J y(if — a)(w — (3)(w 3 + d\w + do) 



where 



« = (8ci)*(^- = (8ci)*(^- - 1); (179) 

OCi DCi 



/ «i bl 3 _ 3 b| Mi (fri+c 2 ) , 
^4 Cl 12c? 4 j ' d °" 8Cljjl 108c? 24c? 8ci ' 



Denote 

F(w) = (w- a)(w - (3){w 3 + d x w + d Q ). (181) 
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Rewrite it as follows: 

F(w) = w 5 + pw 3 + qui 2 + rw + s, (182) 



where 

p = a/3 — G?i(a + 0), q = do — d\ (a + /3), r = a/3d\ — do(a + (3), s = c?oa/3.(183) 



But because of special form of F(w) (see Ea. (|181[l '). we know that the case 4, 
case 5, case 7 and case 12 in section 4 are impossible to occur in the case of 
triple Sine-Gordon equation. Other cases are possible. Thus we can give all 
traveling solutions to triple Sine-Gordon equation according to the cases 1-3, 
case 6 and cases 8-11 in section 4. We omit them for briefly. 

Example 2. All traveling wave solutions to D + 1 dimensional Klein- Gordon 
equation with fourth order nonlinear term 

D 4 

^2u XiXl - u tt = ^bjU 3 . (184) 

i=l j=0 

Take traveling wave transformation u — £i = ^2iLi^i x i + Then 
Eq.jnUJl becomes 

5 

{u'f =a + Y, a ^ u ^ ( 185 ) 

i=l 

where ao is integral constant and 

2 '"' (186) 



(z + iXEf^fc 2 -^ 2 ) 



In order to reduce Ea. (|185|) to the form (|132|l . we take transformations as 
follows: 

i 4 — i i 
w = alu+-a 4 a 5 1 ', £ = a 5 5 £i. (187) 

Then Ea. l|185fl becomes 

±(£-£o) = / dW (188) 

J \/ w 5 + pw 3 + qw 2 + rw + s 



where 



16 2 -f , -f 32 3 _| 

p = ya 4 a 5 + a 3 a 5 , 9 = -_a 4 a 5 + (a 2 - 3a 3 )a 5 , 
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256 4 , _ 4 _^g, 48 2 8 -2 -i 

r = 125 a4 ^ 5 ~ &5 ° ' + 25 a3Ct4a5 ° ~~ 5 a2fl4a 5 + °i a 5 " > 

256 , 4 64 , 16 2 _ 2 4 , . 

S= 3125 a4a5 ~ Y25 a3a4 ° 5 +^° 2a 4 a 5 -gaia^s +"o- ( 189 ) 

According to the cases 1-12 in section 4, it is easy to give the classification of 
all traveling wave solutions to Ea. ll84fl . We omit them for briefly. 



5 The solutions of irrational integrals and appli- 
cations 

There are some mathematical physics equations whose reduced ODE's have the 
form of (u') 2 = j^jy, where g{u) and f(u) are two polynomials. So the classifi- 
cations of their all traveling wave solutions are more complex. 

Example 1. Classification of all single traveling wave solutions to Fujimoto- 
Watanabe equation 

Fujimoto-Watanabe eauation(|37|) reads 

u t = u 3 u xxx + 3u 2 u x u xx + 3eu 2 u x , (190) 

under the traveling wave transformation u = it(£), £ = kx + ujt, the correspond- 
ing reduced ODE is as follows: 

*(e-60=/ , 3 + ud :+ + , (191) 

where 03 = — %,ai = — %r 5 a 2 and ao are two arbitrary constants. For sim- 
plicity, we take the change of variable w = a|u, then the Ea. (|191|) becomes 

±ah-to)=f - .. d W tl d TT > (192) 
J V 10 + d 2 w z + d\w + a 

_ 2 _ 1 

where d 2 = — a 3 3 02, di=a 3 3 a\, do = ao. Denote F(w) — w 3 +d,2W 2 + d\w + 
do, its complete discrimination system is A = — 27(^r + do — ^a) 2 — 4(di — -g 2 ) 3 

d 2 

and Di — di — ^ . There are the following four cases to be discussed: 

Case 1: A = 0,L>i < 0, then we have F(w) = (w - a) 2 (w - 0),ot ^ (3. If 
w > (3, the solutions are as follows: 



dz ( /J(i-i (1 ) = 2 V , a| U -/3+^^=arctan^|L^, a > /?; (193) 

VP - a VP - a 
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±af (£-&) = 2 



(3- 



■An 



a| u — (3 — v 'a — (5 



a>(3; (194) 



a| u — (3 + \/ol — (3 



Case 2: A = 0, Di = 0, then we have F(w) = (w — a) 3 , the solution is as 
follows: 

2 l~. I a 
±af(e-eo) = 2Vaf"-/3- i j ( 195 ) 



a| m — /? 

Case 3: A > 0, D\ < 0, then F(w) = (w — a)(w — (3){w — 7), we suppose 
that a < (3 < 7. It is easy to see the corresponding integral can be expressed 
by the second kind of elliptic integrals. 

Case 4: A < 0, then we have F(w) = (w — a)(w 2 + pw + q), p 2 — 4g < 0. 
Then the corresponding integral can be expressed by the second kind of elliptic 
integrals. 

Example 2. The classification of all single traveling wave solutions to cou- 
pled Harry-Dym eauationf [38j) 

v, = -u,v-\„ I + ev,. (197) 



Through traveling wave transformation and integrations, we obtain 

Co 



(198) 



k 2 u + ro 



{\ku + oj — fee) 2 ' 

A 2 4- n, in 4- r„\ 4- n„ 

(199) 



, 2 _ 2 (u + r ) 4 + a 3 (u + r ) 3 + a 2 (u + r ) 2 + ai(n + r ) + a 

( M J — 7.9 ... , „ 



where r = ^-2e,a 3 = -4(^-e,a = ^ci. Denote F(n) = (M + r ) 4 + a 3 (u + 
ro) 3 + a2(u + ro) 2 + a%(u + J"o) + aoj according to the section 4, the complete 
discrimination system for the fourth order polynomial F(u) is as follows: 



D x = 4, D 2 = -p, D 3 = 8rp - 2p 3 - 9q 2 , 
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F 2 = 9p 2 - 32pr, (200) 
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L» 4 = 4p 4 r - p 3 q 2 + 36prq 2 ~ 32r 2 p 2 - —q A + 64r 



2G 



where p = a 2 , q = - ^ + a i); r = a o - ^ + ^ - |§- There are the 
following cases to be discussed. 

Case 1:D 4 = 0, D 3 = 0, D 2 < 0. Then we have 

F{u) = (u + r ) 2 + l(u + r ) + s 2 , (201) 
where I, s are real numbers, I 2 — 4s 2 < 0. We have 



± (£ _ £ ) — 1 ln " + r oT vgs - + r + s 
2v2s— T ?i + ro ± \/2s — + r o + s 

+ -^{arctan V ^£^T + ^ A^T7jTV^ l y 



Case 2: L> 4 = 0, L> 3 = 0, D 2 = 0. Then we have 

F(u) = (u + r a - a) 4 . (203) 

If a > 0,we have 

If a < 0,we have 

j_ ic c \ 1 +. \ /u + r o \/u + r a 

± — 4q) = arctan — == . (205) 

v/— a v— a u + ro — a 



Case 3: D 4 = 0, D 3 = 0, D 2 >0,E 2 = 0. Then we have 
F(u) = (u + r - af(u + r - (3) 2 , 
where a, (3 are real numbers, a > (3. If a > (3 > 0, we have 

.y/u + r - v 7 ^. 



± (a + /?)(£-&) 



/a ln 



V^lnj 



Vu + ro + Vfi 



If > a > j3, we have 

± (a + /?)(£- £o) = 2 arctan 
If a > > (3, we have 



±(a + /?)(£-£o) = \/^m 



V" + ni + \/a 



— 2y/—(3 arctan 



— 2\J—(3 arctan 



(206) 
(207) 

(208) 

(209) 



Case 4. In all other cases such as F(u) = (u — a) 2 (u — j3)(u — 7) and so on, 
the corresponding solutions can be expressed with elliptic functions and hyper- 
elliptic functions. We omit them for briefly. 
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6 Classifications of solutions to some more gen- 
eral nonlinear differential equations 

There are many nonlinear differential equations which can be reduced to an 
integrable ODE as follows 

u"(Z)+p(u)(u'(Z)) 2 + q(u) = 0. (210) 

These equations include High order KdV equation( 50 ) 

u t + u x + auu x + {3u xxx + a 2 piu 2 u x + a(3(p 2 uu xxx + p 3 u x u xx ) = 0; (211) 



Coupled KdV equation(jSJ) 

u t + eiv x + e 2 u 2 u x + e 3 u xxx = 0, (212) 
v t + r)i(uv) x + i] 2 vv x = 0; (213) 

Coupled mKdV equation([S2j) 

u t + Eivv x + e 2 uu x + e 3 u xxx = 0, (214) 
v t + i]i (uv) x + rj 2 vv x = 0, (215) 

and so on. We give the general solution of the Ea. (|210|) in the following: 

Lemma 53]: The general solution of ODE 121UI) is as follows: 

f du , 

±(£ - &) = / , (216) 
J wexp(— 2 J p(u)du)[c — 2 J q(u) exp(2 J p(u)du)du] 

where c and £o are t w0 arbitrary constants. 

We take the coupled KdV equation as an example to illustrate our method 
for shortly. Other equations can be dealt with similarly. 

Example: Classification of all single traveling wave solutions to coupled 
KdV equation. 

Take the traveling wave transformation as« = «((),j) = i)((), £ = kx + ut, 
its reduced equations are as follows: 

7 /c V + ^-u 2 - ^-v 2 + luu + c = 0, (217) 
where cq and c\ are two integral constants. There are two cases to be discussed. 
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Case 1. c\ = 0. We have u = — — and w satisfies the following 
equation 

{v'f = b 3 v 3 + b 2 v 2 + b x v + 6 , (219) 

where 

2p ak 0kS 2 _ 06m -Sup 2p 0co 2 2^ 2 

b3 -^ { T + Tv } ' h ~ 7 fcV ' bl -W6 { Jp T ~Tp~ + oh 

(220) 

and bo is an arbitrary constant. Denote e?2 = ttj^i = TT'^ = j 2 - F(v) = 

w^^ + dxw + do, and A = -27d§-4(-^) 3 (-%)-\.Di = - ^)~*. 
Here A and Z?i make up a complete discrimination system for F(v).We give the 
classification of all solutions to the Ea. H219fl as follows: 

Case 1.1: If A = 0, D x < 0, then we have F(v) = (v - a) 2 (v - 0),a^ 0. If 
v > 0, the solutions are as follows: 







V 2 




f b 3 ^/a - 





v 2 


,hV0- 


01 , 



(£-£o))+/8, a>0; (221) 

(£-&)) + A «>/3; (222) 
(£-&)) + <*, a</3. (223) 



Case 1.2: If A = 0, Z?i = 0, then we have -F(w) = (w — a) 3 , the solution is 
as follows: 

v = 4b 3 (H - to)- 2 + a. (224) 



Case 1.3: If A > 0, D\ < 0, then F(v) = (v — a)(v — 0)(v — 7), we suppose 
that a < < 7. When a < v < 0, we have 

v = a + (0 - a) sn 2 ( 53V ^" (g - fr)),m); (225) 
when u > 7, we have 



7 _ ^ sn 2 ( j3v5E H ( ^_ go )) ;m) 



w = ' ' . X2 " w/ " (226) 
cn 2 (^^(e-^o)),m) 



where 



,2 _ B-a 
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Case 1.4: If A < 0, then F(v) = (v - a)(v 2 +pv + q), p 2 -Aq < 0. We have 



v = a-yja z +pa + q^ ^ J, (227) 

1 + cn((a 2 + pa + q) J ((£ — £ )j m) 

where m 2 = ^(1 , " +? ). 

Case 2. ci ^ 0. According to the lemma, then u satisfies the following 
integral form: 

± ({- & )-/ , f--' d ° . (228, 
J vyagv 8 -I h aiv + ao 

where ai,-- - , as can be expressed by the corresponding parameters. Denote 
F(v) = agv s + ■ • • + ciiv + ao, its complete discrimination system is too complex 
to write. We omit these solutions for briefly. 



7 Discussions and conclusions 

In Ref. JHjj Fan proposed an expanse method to seek for the single traveling 
wave solutions to nonlinear partial differential equations. His key idea is to 
assume that the solution has a polynomial form u = a n <p n + ■ + aic/) + ag of a 
function <j> which satisfies an ODE as follows: 

(j)' = ±Vc r (f> r + • + ci0 + c o , (229) 

According to the balance principle, a relation of n and r can be obtained. In 
his paper, Fan take a special case of r = 4 as an example to apply his method. 
In the end of his paper, he said "we only have investigated a special case when 
r = 4. the proposed method can be extended to the cases when r > 4. The 
details for those cases will be investigated in our future works" . In fact, there 
are two problems need to answer. One problem is whether or no there exist new 
solutions in the cases of r > 4 for some equations. The second problem is how 
to prove some solutions are new in the cases of r > 4 if the equation considered 
can't be directly reduced to an integral form. 

According to the above analysis to direct integral method ,we can make a 
conclusion that there aren't new solutions in the cases of r > 4 for those equa- 
tions which can be directly reduced to an integral form. Thus we have answered 
the Fan's first problem. 



Direct integral method is very natural and simple for solving traveling wave 
solutions to nonlinear partial differential equations. But it is just complete 
discrimination system for polynomial make direct integral method becomes an 
efficient and powerful method. A lot of nonlinear mathematical physics equa- 
tions can be dealt with by direct integral method, their all traveling wave so- 
lutions can be classified. These complete results can't be obtained by using of 
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other methods. On the other hand, a lot of nonlinear evolution equations such 
as (l+l)-dimcnsional Boussinesq equation and (3+l)-dimensional KP equation 
and so on, can be reduced to integral forms by setting some integral constants 
to zeroes. For those equations, we can not give their all traveling wave solutions, 
but we can still obtain their abundant exact solutions. For example, we consider 
the (2+l)-dimensional boussinesq equation 

Utt U X x ~\~ 3(li *)xx Hxxxx ^yy — 0- (230) 

Instituting traveling wave transformation u = = k\x + k 2 y + uit into 

Ea. l|23U[l and integrating twice yield 

k\u = {uj 2 -k\- k\)u + 3klu 2 + c 2 £ + c x . (231) 

where c\ and C2 are two arbitrary constants. If we take ci = 0, then we can 
apply the direct integral method to it. In fact, integrating Ea. (|231|) once yields 

±*i(£ " Co) - / = ( 2 , 2 d " ;2N = (232) 

where cq is an arbitrary constant. According to Ref.[19], we can easily to give 
classification of all solutions to Ea. l|232|l . Thus it is easy to obtain abundant 
traveling wave solutions to Eq. (|230ll . For briefly, we omit them. 
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